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MAGNETIC CONTROL OF THE FLOW OF HOT COMBUSTION GASES IN 
HYDROGEN-OXYGEN COMBUSTION AT 14 ATMOSPHERES 
by Ronald Razner 
Lewis Research Center 

SUMMARY 

The feasibility of using electromagnetic forces for the purpose of directing the flow 
of hot combustion gases was investigated in the combustion of gaseous hydrogen and liquid 
oxygen at a pressure of 14 atmospheres. 

Radiation between 2500 and 3500 angstroms emitted in the combustion process was 
used to detect changes in the extent of spreading of a single jet of liquid oxygen reacting 
in an atmosphere of gaseous hydrogen. This radiation is known to consist of hydroxyl 
radical and oxygen bands plus a continuum. 

In general, the exertion of electromagnetic forces of a useful magnitude depended on 
the existence of a suitably high ion population in the combustion gases. When this popula- 
tion was sufficiently large, the velocity of the gases coupled with an externally applied 
magnetic field of 0. 0200 to 0. 0380 weber per square meter to exert an influence of prac- 
tical magnitude on the directional flow of the gases. 

Calculated estimates are made and discussed to compare the effectiveness of two 
methods for increasing the ion population in hot gases of a type encountered in hydrogen- 
oxygen combustion. 


INTRODUCTION 

A persistent problem in the design of chemical combustors involves the control of 
damage to mechanical surfaces resulting from contact with hot combustion gases. This 
problem increases in complexity with increasing temperature and pressure. Conse- 
quently, as a part of effective combustor design, investigating all possible methods of 
heat-transfer control is desirable. 

To date the most serious heat-transfer problems known occur in a field quite far re- 
moved from chemical combustion. These problems are in areas of plasma physics con- 



cerned with controlled thermonuclear fusion, where it is 
suggested that heat transfer can be controlled by means of 
magnetic confinement. In recent years extensive theoreti- 
cal and experimental research has been done in this area. 

While the physical conditions of high-pressure chemi- 
cal combustion are quite different from those encountered 
in the physics of highly ionized and fully ionized plasmas, a 
consideration of the electrical properties of hot combustion 
gases (ref. 1) suggests that certain aspects of the principle 
of magnetic confinement might still be applied. The practi- 
cability of this concept would be established if it could be 
shown that directional flows can be significantly influenced 
by applied magnetic fields of reasonable intensity. 

The basic problem is depicted schematically in fig- 
ure 1. A volume element of combustible material Q near a mechanical surface is sub- 
ject to a sudden expansion through the process of combustion. It is shown in appendix A 
that if the expanding gas has a sufficiently high electric conductivity, the presence of a 
magnetic field B parallel to the surface will set up forces that result in a preferred di- 
rection of gas expansion parallel to the surface. This effect increases as the electric 
conductivity increases. 

The effect of an externally applied magnetic field on the expansion of a single jet of 
liquid oxygen burning in an atmosphere of gaseous hydrogen was studied at a pressure of 
14 atmospheres in a rocket combustor. A provision was made for increasing the ion 
population in the combustor by using the moderate (initial) electric conductivity of the 
burning jet to support a small electric current discharge. 

Spectral radiation of known origin (ref. 2) emitted in hydrogen- oxygen combustion 
was used photographically to detect and to study the desired effects. Since the radiation 
intensity increases as a function of temperature, any influence on the flow stream de- 
tected by using this radiation can be interpreted as an influence on the high- temperature 
portions of the gases. 

An estimate was made of the lower limit of the ion population produced by an electric 
current flowing in the combustion gases. This limit is compared with calculated esti- 
mates of thermal ion populations produced by the technique of seeding. 



Figure 1. - Visualization of volume element 
of fluid in a flow stream with impressed 
magnetic field. 


APPARATUS AND PROCEDURE 

The combustor (fig. 2) was constructed in sections consisting of an injector, boron 
nitride electrical insulators, magnetic coils, window section, and nozzle. All parts were 
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made of copper, except for electrical insulators. Nylon pipe connectors insulated the 
combustor electrically from all external liquid- and gas- supply lines. Silicon rubber O- 
rings provided pressure seals between adjacent combustor sections. 

The injector consisted of a single jet 0. 100 inch in diameter alined to inject a stream 
of liquid oxygen down the cylindrical axis of the combustor. This jet was surrounded by 
a parallel injected flow of gaseous hydrogen. A 600- volt electric potential placed between 
the injector and the nozzle provided an electric current discharge along the burning jet 
for the purpose of increasing the ion population in the hot gases. The injector was made 
the cathode and the nozzle, the anode. Current flowing in the jet was monitored by re- 
cording the voltage drop across a 1-ohm series resistor. 

Two magnetic coil windings, each consisting of 100 turns of number 14 copper mag- 
net wire, produced a magnetic field parallel to the axis of the combustor. These were 
wound on 4-inch outside diameter copper engine sections. The width of each winding was 
5/8 inch, and the coil ends were separated on the combustor by 5^ inches. Each winding 
was powered by 30. 0-amperes direct current supplied from two 12-volt storage batteries 
connected in series. This provided a (measured) magnetic field on the chamber axis which 
varied from 0. 0200 weber per square meter halfway between the coils to 0. 0380±0. 0005 
weber per square meter in the plane of the coils. 

A 2- by 1/4- inch viewing port equipped with quartz windows provided a cross- 
sectional view of the 2 -inch-diameter chamber at right angles to the chamber axis. This 
port was placed halfway between the magnet coils where the magnetic field had an inten- 
sity of 0. 0200 weber per square meter. A helium purge was used to keep the recessed 
ports and windows free from obstruction by absorbing gases and solid particles resulting 
from erosion of combustor surfaces. 

Radiation from the viewing port was focused on the slit of a grating spectrograph 
operating in the spectral range of 2500 to 4500 angstroms in the first order. The 2-inch 
length of the port was focused parallel to the slit of the spectrograph so that long spectral 
lines were produced (fig. 3) with points along their lengths corresponding to points on the 
viewing cross section of the combustion chamber. The slit width was set at 30 microns. 
This technique provided a means of simultaneously studying effects on the spreading of 
the burning oxygen jet and the emission spectra associated with the combustion. Films 
were processed and spectral intensity data obtained by methods discussed in reference 2. 

The combustor was operated at an oxidant- fuel-weight ratio (o/f) of 3. 88 with a 
maximum deviation from this value of 0. 07 for any single run. The average total weight 
flow for all runs was 0. 549 pound per second with a maximum deviation for any single 
run of 0. 007 pound per second. A chamber pressure of 14. 0 atmospheres was maintained 
to within 0. 1 atmosphere on all runs. All electrical sources were turned on and the 
spectrograph slit opened before the start of the combustion, which was 3 seconds in dura- 
tion. 
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(a) Combustion in absence of magnetic field and electric 
current flow. 



(b) Combustion in presence of impressed magnetic field 
but without electric current flow. 



(c) Combustion in absence of impressed magnetic field 
but with electric current flowing through the gases. 



id ) Combustion in presence of both impressed magnetic 
field and electric current flow in the gases. 

Figure 3, - Hydroxyl radical spectra showing extent of spreading at viewing 
position of burning oxygen jet (band head at 3064 Ah 



With these conditions four types of 
experiments were carried out (herein- 
after referred to as experiments a, b, c, 
and d, respectively). These were (a) 
control runs consisting of the combustion 
in the absence of impressed fields and 
electric current flow; (b) combustion in 
the presence of a magnetic field but with- 
out electric current flow in the burning 
jet; (c) combustion in the absence of a 
magnetic field but with electric current 
flow produced by the 6 00- volt applied po- 
tential; (d) combustion in the presence 
of both a magnetic field and electric cur- 
rent flow. 

EXPERIMENTAL RESULTS 
Spectra 

A typical set of spectra correspond- 
ing to experiments (a) to (d) is shown in 
figure 3. The predominant feature is the 
molecular band head of the (O, O) branch 

2 2 -T-T- 

in the 53 — JJ~ electronic transition 
of the OH radical, which occurs at 3064 
angstroms. Other vibrational branches 
for this transition, including (1,0) and 
(0, 1), were present on the films but 
were too weak to show up in the prints. 
The qualitative nature of the radiation 
remained the same throughout all the ex- 
periments. In general, effects of self- 
absorption, broadening, background, and 
overlapping spectra of other molecular 
bands besides those due to OH were ob- 
tained as in previous radiation studies 
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(ref. 2). Difficulties presented by these interferences were such that no attempt was 
made to interpret the details of rotational fine structure. A comparison of the general 
intensity levels of the (O, O), (1,0), and (0, 1) vibrational branches in OH indicated that 
the (0, 1) intensity may have increased relative to the others whenever an electric current 
was flowing in the gases. (Within the limits of experimental error the (1, O) and (O, O) in- 
tensity levels, as photographically recorded, remained the same from film to film. ) How- 
ever, this effect was difficult to distinguish because of the many errors involved in com- 
paring photographic spectral intensities over the required wavelength intervals. Conse- 
quently, while relative vibrational intensity changes in the OH radiation would not be sur- 
prising, little confidence or importance was placed on such changes in this work, and 
their investigation was considered outside the scope of the present report. 


Electrical Properties of Burning Jet 

Whenever the 600- volt potential was placed between the injector and the nozzle, a 
current flow that varied erratically between 30 and 200 milliamperes was obtained during 
the entire duration of combustion. The existence of actual magnetic forces on the flow 
stream (as is evident from the data given in the section Electromagnetic Effects on Flow 
Stream) established that the current was flowing through the combustion gases and could 
not be attributed to a leakage path by way of the copper walls of the combustor. On the 
basis of Ohm’s law, the observed currents corresponded to an overall electrical resis- 
tance of the burning jet of 3000 to 20 000 ohms. The lower value corresponded to rather 
peaked current surges during the combustion and indicated tendencies toward momentary 
high degrees of electrical breakdown. 


Electromagnetic Effects on Flow Stream 

The extent of spreading of the combustion zones in experiments (a) to (d) (as indi- 
cated by the emitted radiation) is shown in figures 3 and 4. In figure 4 the fraction of 
maximum intensity l/l^ ov of the (O, O) band head of OH for each spectrum of figure 3 
was plotted against position on the chamber cross section. Figures 4(b) and (c) show a 
slight narrowing when a magnetic field or electric current was applied separately. How- 
ever, figure 4(d) shows that a marked narrowing occurred when a combination of these 
two influences was applied. In this case the combustion zone was contracted to about 
one-half of that obtained with no external influences. (A few runs were made where the 
discharge did not exceed 150 mA because of film deposition on the nozzle. In such in- 
stances the magnitude of the constriction was somewhat reduced. ) 
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DISCUSSION 


The results shown in figure 4 indicate 
that magnetic control of the directional 
flow of combustion gases in hydrogen- 
oxygen combustion at 14. 0 atmospheres 
pressure can be achieved provided that a 
sufficiently high ion population and electri- 
cal conductivity can be maintained. It ap- 
pears that in the combustion system with 
only the magnetic field applied these con- 
ditions might not be realized. 

Two simple techniques that can be 
used for increasing the ion population in a 
hot gas are seeding and the application of 
an electric discharge. Seeding increases 
the thermal ion population by adding atomic 
and molecular species with relatively low 
ionization potentials to the system. If this 
technique is used, one must choose a seed- 
ing material with a low ionization potential, 
and also one that will not react chemically 
to form molecular species with high disso- 
ciation energies. Such occurrence would 
defeat the aim of reducing the ionization 
potential of the system. Consequently, 
the method of seeding requires a consider- 
ation of the chemistry of each individual combustion system. Table I is a summary of 
thermal ion populations calculated for three ionization potentials by assuming the validity 
of the Saha equation (ref. 3, p. 92) and the ideal gas law. From this table it is seen how 
the ion population changes with ionization potential, temperature, and pressure, and that 
the above conclusion about the ion population in the combustion system, taken by itself, 
might be expected for the temperature range of interest (3000° to 4000° K). 

In this work a need for any detailed chemical considerations has been avoided by 
passing an electric discharge through the gases. This technique increased the ion popu- 
lation and electrical conductivity by producing and sustaining a state of electric break- 
down in the gases. The ion population and conductivity produced were sufficient to per- 
mit a magnetic deflection of the flow stream. In particular it should be noted that the 
electrical currents and magnetic field strengths needed to produce this effect were not 



(d) Intensity corresponding to figure 3(d). 

Figure 4. - Intensity of ultraviolet OH emission across viewing 
cross section. 
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TABLE I. - TOTAL THERMAL ION POPULATIONS FOR SINGLY CHARGED IONS 

[M — M + + e'. ] 


Temper- 



Ion population 

a N , ions/m ^ 



ature, 









°K 




Pressure, 

P, atm 





1 

10 

20 

30 

40 

50 

60 

70 


Ionization potential for hydro gen -oxygen combustion, ^ v., 12.5 


1 000 

1.3X10' 6 

4. lxlO' 6 

5. 9X10' 6 

7.0X10' 6 

8.2X10* 6 

8.8X10' 6 

9. 7X10' 6 

lOxlO -6 

2 000 

0. 93X10 10 

2.9X10 10 

4. lxl 0 10 

5. lxlO 40 

5. 9X10 40 

6.6X10 10 

7. OxlO 40 

7.7X10 40 

3 000 

1.8X10 15 

5. 9x1 0 4 5 

8.3X10 15 

lOxlO 15 

llxlO 45 

13X10 45 

14X1 0 4 5 

15X10 45 

4 000 

0. 84x1 0 18 

2.6X10 18 

3.7X10 18 

4.6X10 48 

5. 3xl0 48 

5. 9X1 0 48 

6. 4xl0 18 

6. 9X10 48 

5 000 

0. 35X10 20 

1.1X10 20 

1. 5x1 0 20 

1. 9X1 0 20 

2. 2x1 0 20 

2. 5xl0 20 

2. 6xl0 20 

2. 9xl 0 20 

6 000 

0. 39x1 0 21 

1. 3x1 0 21 

1. 8x1 0 21 

2. 2x1 0 24 

2. 5xl 0 24 

2. 8x1 0 21 

3. lxl 0 24 

3. 4x1 0 24 

7 000 

0. 23X10 22 

0. 73X10 22 

l.OxlO 22 

1. 2x1 0 22 

1.4X10 22 

1. 6x1 0 22 

1. 8xl 0 22 

1. 9xl 0 22 

8 000 

0. 88X10 22 

2. 8x1 0 22 

4. 0x1 0 22 

4. 8xl0 22 

5. 6xl0 22 

6. 2xl0 22 

6.8X10 22 

7. 3X1 0 22 

9 000 

0. 24x1 0 23 

0. 80xl 0 23 

1. lxlO 23 

1. 4x1 0 23 

1. 6x1 0 23 

1. 8x1 0 23 

2. OxlO 23 

2. OxlO 23 

10 000 

0. 59x1 0 23 

1. 9X10 23 

2. 6x1 0 23 

3. 2x1 0 23 

3.7X10 23 

4. lxlO 23 

5. OxlO 23 

4. 9xl0 23 



Intermediate ionization potential, v-, 
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1 000 

0. 32X10 11 

1.0X10 11 

1. 5x1 0 4 4 

1.8X10 44 

2.0X10 44 

2.3X10 44 

2. 5x1 0 11 

2. 8x1 0 11 

2 000 

1. 5xl0 18 

4. 6X10 18 

6. 6X10 18 

7. 9X10 48 

9.4X10 48 

lOxlO 18 

llxlO 48 

12X10 48 

3 000 

0. 54xl0 24 

1.7X10 21 

2.4X10 24 

3. lxlO 21 

3. 5xl0 21 

3. 9X10 24 

4. lxlO 24 

4. 4xl0 21 

4 000 

1. 1X10 22 

3.4X10 22 

4. 8X10 22 

5.8X10 22 

6.8X10 22 

7. 5xl0 22 

8. lxlO 22 

9. OxlO 22 

5 000 

0.64X10 23 

2.0X10 23 

2. 9xl0 23 

3. 5xl0 23 

4. 0X10 23 

4. 6xl0 23 

4. 9xl0 23 

5. 3xl0 23 

6 000 

0.23X10 24 

0.69X10 24 

0. 98X10 24 

1.2X10 24 

1.4X10 24 

1.6X10 24 

1. 8xl0 24 

1. 9X10 24 

7 000 

0. 50xl0 24 

1. 7xl0 24 

2.3X10 24 

2.8X10 24 

3.3X10 24 

3.7X10 24 

4. OxlO 24 

4. 2xl0 24 

8 000 

0. 99xl0 24 

3. 1X10 24 

4.4X10 24 

5.4X10 24 

6.2X10 24 

7. OxlO 24 

7. 6xl0 24 

8. 2xl0 24 

9 000 

1.2X10 24 

4. 9X10 24 

7.2X10 24 

8.8X10 24 

lOxlO 24 

llxlO 24 

13xl0 24 

14X10 24 

10 000 

0. 13X10 25 

0.70X10 25 

l.OxlO 25 

1.3X10 25 

1. 5xl0 25 

1.8X10 25 

1. 9xl0 25 

2.2X10 25 



Ionization potential for pure cesium gas, v., 3. 87 



1 000 

0.79X10 16 

2. 5xl0 16 

3. 5xl0 46 

4. 3xl0 16 

5. OxlO 16 

5. 6X10 46 

6. lxlO 46 

6.6X10 46 

2 000 

0.72X10 21 

2.3X10 21 

3.2xl0 24 

4.0X10 24 

4. 4xl0 21 

5. lxlO 24 

5.7X10 24 

6.2X10 24 

3 000 

0. 34X10 23 

1. lxlO 23 

1.6X10 23 

1. 9X10 23 

2.2X10 23 

2.4X10 23 

2. 6xl0 23 

2.8X10 23 

4 000 

0.24X10 24 

0. 73X10 24 

l.OxlO 24 

1.3X10 24 

1. 5xl0 24 

1. 6X10 24 

1. 8xl0 24 

2. OxlO 24 

5 000 

0.76X10 24 

2.4X10 24 

3. 5xl0 24 

4.2X10 24 

4. 8X10 24 

5. 4xl0 24 

5. 8xl0 24 

6.6X10 24 

6 000 

1.8X10 24 

5.4X10 24 

7. 8xl0 24 

9. 5xl0 24 

llxlO 24 

12X10 24 

13X10 24 

14X10 24 

7 000 

0. 17X10 25 

0. 86X10 25 

1.2X10 25 

1.6X10 25 

1.8X10 25 

2. lxlO 25 

2. 3xl0 25 

2. 5xl0 25 

8 000 

0. 17X10 25 

1. lxlO 25 

1.8X10 25 

2.3X10 25 

2.7X10 25 

3. OxlO 25 

3.4X10 25 

3.8X10 25 

9 000 

0. 16X10 25 

1.3X10 25 

2.2X10 25 

2.9X10 25 

3.4X10 25 

4. OxlO 25 

4. 5xl0 25 

4. 9xl0 25 

10 000 

0. 14X10 25 

1.3X10 25 

2.3X10 25 

3.2X10 25 

4. OxlO 25 

4.7X10 25 

5. 4xl0 25 

6. OxlO 25 


Ion population is N Q = 2^ where log[x 2 /(l - x 2 )jp = (-5050 Vj/T) +2.5 log T - 6. 5. 

^Estimated from table 4. 4, ref. 3, p. 83. 



unduly high. On the basis of total weight flow, the electric current flow in the gases 
amounted to 0. 055 to 0. 36 coulomb per pound, or in terms of power, 33 to 220 watt- 
seconds per pound. The influence of these currents can be seen by considering the fluid 
system in terms of its classification as a plasma. According to accepted usage (ref. 4, 
p. 272), a fluid can be classified as a plasma if L » h, where L is the characteristic 
length of a volume element of fluid Q and h is the Debye shielding length given by 


h = 




(meters) 


The last expression gives h in meters if T is in degrees Kelvin and N Q is in ions per 
cubic meter. (All symbols are defined in appendix B. ) For our purposes, this quantity 
serves as a criterion distinguishing macroscopic systems (which can be treated by meth- 
ods of macroscopic fluid mechanics) from those requiring a more detailed "microscopic" 
treatment. From table I, it is seen that if one relied only on thermal ionization to deter- 
mine the ion population in hydrogen- oxygen combustion, N would be about 7x10*® ions 
per cubic meter at 3000 K and about 3x10 ions per cubic meter at 4000° K (for 14. 0- 
atm pressure), which corresponds to a range in h of 3x10”® to 2x10”® meter. On the 

other hand, it is calculated in appendix A that the electric discharge maintained a mini- 

1 7 

mum ion population of 3x10 ions per cubic meter (assuming a burning temperature 
range of 3000° to 3500° IQ, which corresponds to a value for h of 5x10"® meter or less 
during the combustion. 

The present combustor achieved a rather low c* efficiency (ref. 5) of around 50 per- 
cent, with a slight lowering (about 1^ percent) observed in experiments of type d. The 
lowering of c* efficiency was probably due to a reduction in the mixing efficiency in the 
presence of magnetic confinement. While this result is not surprising for the combustor 
configuration used, it is of little consequence for this study. The injector system was 
not chosen for the achievement of high c* efficiencies but rather as an ideal observa- 
tional system in which magnetic influences on a burning flow system could be studied as 
depicted in figure 1 (p. 2). 

In addition to modifying heat transfer, magnetic forces could likely be used to con- 
trol the amplitude of hydrodynamic flow oscillations in a combustor by exerting direc- 
tional damping influences on the flow system. 


SUMMARY OF RESULTS 

A study of the influence of a magnetic field on the directional flow of hot gases in 
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hydrogen- oxygen combustion at a pressure of 14. 0 atmospheres yielded the following re- 
sults 

1. An applied magnetic field of 0. 02 weber per square meter was sufficient to influ- 
ence flow direction in a combustion flow stream at several atmospheres pressure pro- 
vided that the gases had a sufficiently high electrical conductivity. 

2. Apparently the electrical conductivity required for this influence is higher than 
that normally existing in hydrogen- oxygen combustion taking place in the temperature 
range of 3000° to 3500° K at 14. 0 atmospheres pressure. However, this can be in- 
creased sufficiently by allowing a small electric current discharge to flow in the combus- 
tion system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 17, 1965. 


10 



APPENDIX A 


ELECTROMAGNETIC PROPERTIES OF COMBUSTION FLUID SYSTEM 

A theoretical consideration of electromagnetic forces that can act on a fluid system 
involves a study of the equations of magnetic fluid dynamics. For a macroscopic formu- 
lation (which will be used) these equations are obtained by combining Maxwell's equations 
with the equations of fluid dynamics and are discussed in some detail elsewhere (ref. 6, 
ch. 6; ref. 7, ch. VIII; ref. 4, ch. 14). Other possibilities include microscopic formu- 
lations using the principles of kinetic theory. In general, these may be quite complicated, 
and the physical conditions to be introduced will vary greatly with types of electrical prop- 
erties, as well as with the nature of fluid motion. There are several factors in the case 
presented herein which seriously restrict our ability to attempt any detailed theoretical 
treatment (even in the macroscopic sense). Not the least of these is a lack of quantitative 
understanding of any but bulk electrical properties of combustion gases undergoing chemi- 
cal reaction; this lack in turn is coupled with problems of mixing and nonsteady flows in 
the combustor. Consequently, in this appendix we shall proceed mostly by analogy to an 
idealized problem that clearly exhibits the physical conditions necessary to achieve mag- 
netic control of flow direction and make certain general observations about the nature of 
this fluid system. 

Consider a volume element of fluid Q with characteristic length L, electrical con- 
ductivity a, and fluid density p, moving with velocity u in the presence of local electric 
and magnetic fields E and B, respectively. (The magnetic permeability ju will be 
taken as 1, thus treating any paramagnetic properties such as those of molecular oxygen 
as of secondary importance. ) The fluid can be considered a conductor if the characteris- 
tic electrical relaxation time t (ref. 4, p. 140) is small compared with characteristic 
times of flow fluctuations. (This will be sufficiently good for our purposes if the fre- 
quency defined by u> = 1/t is greater than about 1 Mc/sec.) 

In general, the motion of a conductor in a magnetic field will cause electric currents 
to flow, resulting in complex interactions with all electric and magnetic fields that may 
be present and with the flow itself. The total electric current density j in volume ele- 
ment Q is related to the fields E and B through Ohm's law: 

J = ct(E + u x B) (Bl) 

In addition, the magnetic field is subject to a divergence condition 

V • B = 0 (B2) 
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and a time rate of change governed by Faraday’s law 


V x E = - — (B3) 

at' 


while the general fluid motion itself conforms to the equation of continuity 


& + V • (pu ) = 0 (B4) 

at 


Now consider the condition where a is sufficiently large (see following discussion) that 


1-0 (B5) 

cr 


Then 


E = -u x B or V X E = -V X (u X B) (B6) 

Substituting equation (6) into (3), expanding the result by vector identity, and noting condi- 
tion (2) give 

— = (B • V)u - (u ■ V)B - BV • u (B7) 

at 

Using the condition V • u = -(l/p)[(3p/at) - u • Vp] from equation (4) and noting that the 
operator a/at + u • V h D/Dt is the "substantial derivative" allows equation (7) to be 
written as 


£(I) = (I)<b.v ) u - (BS) 

On the other hand, it is easily shown (ref. 7, p. 216) that if £ represents an element of 
length on a fluid line, its rate of change is given by 

— (f) = (T • V)u (B9) 

Dt 

Consequently, when the conductivity conditions for the derivation of equation (8) hold, 
both B/p and jf experience rates of change governed by similar formulas. If, as is 
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presently the case, an external field B is imposed on a reaction chamber prior to the 
presence of fluids in motion, then to the extent that B can remain unaltered upon injec- 
tion of a conducting fluid, a change in L must correspond to a change in p in the direc- 
tion of B. 

A point must now be made concerning the application of this derivation that requires 
a more quantitative view of the fluid system. For this purpose three quantities must be 
considered: the Reynolds number Re, the magnetic interaction parameter N, and the 
magnetic Reynolds number Re m (ref. 6, p. 172). These are given by 

Re = — (BIO) 

v 


where v represents the kinematic viscosity, 


N = 


ctB 2 L 

pu 


and 


Re m = 


(Bll) 


(B12) 


The product of equations (10) and (11) defines another quantity, the Hartmann number, 
given by 


H 2 - B 2 L 2 g 
P v 


(B13) 


and represents a measure of magnetic forces as compared with viscous forces. The 
magnetic interaction parameter represents the ratio of u x B forces to inertial forces, 
whereas the magnetic Reynolds number measures the distortion that the fluid motion pro- 
duces on the magnetic field. For magnetic forces to dominate viscous and inertial forces, 
H > 1 and N > 1, whereas for the magnetic field to remain unaltered by the flow 

Re m « *■ 

We will now attempt to establish maximum ranges for Re, N, Re m , and H by taking 
fluid dimensions to be characterized either by the full chamber length (for the assumption 
of a completely continuous fluid) or by the dimensions of a small pocket of combustion 
gases (about 10”^ m), with the electrical conditions as deduced from the overall electri- 
cal resistance of the combustion system. The pertinent quantities are then obtained as 
follows: 
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Viscosity (estimated from ref. 8), kg/sec-m 7xl0“ 3 to 8xl0“ 3 

Average molecular weight (estimated for equilibrium composition at 

o/f=3.88), kg/mole lXlO -3 

Average molar population (from ideal gas law at 14 atm and 3000° to 

3500° K), moles/m 3 49. 0 to 57. 0 

3 

Average density, kg/m 0. 49 to 0. 57 

Kinematic viscosity (viscosity/density), v, m 3 /sec 1. 2x10“^ to 1. 6x10“^ 

Average weight flow, kg/sec 0.25 

2 -3-3 

Area (1/2 chamber radius to full chamber radius), m 2x10 to 8x10 

— 4 

Length (from gas pocket size to full chamber length), m . . . 10 to 0. 28 

Average velocity (free stream), weight flow/(density)(area), m/sec 54 to 250 

Conductivity, length/(resistance)(area), (for length, 

0.28 m), mho/m 1.7x10 3 to 4.7X10 3 


From these quantities maximum ranges for Re, N, Re m , and H are obtained: 

3.4X10 1 <Re< 5. 8X10 5 
4. 9xl0' 13 < N < 2. OxlO -7 
9. 2xl0 -6 < Re m < 3. 3x10° 

4. lxlO -6 <H < 3. 5X10“ 1 

On the basis of these calculations both the Hartmann number and the magnetic inter- 
action parameter are appreciably less than unity, and hence the magnetic forces cannot 
dominate the viscous and inertial forces present in the combustor. It is evident that the 
macroscopic considerations used herein are not adequate for explaining the observed flow 
constriction caused by the current flow and the applied magnetic field. Valid conclusions 
cannot be obtained, therefore, from calculations based on free- stream flow conditions 
and the assumption of fluid continuity. This is particularly the case if electric conductiv- 
ity is estimated from the electric current drawn by the overall flow stream. In this re- 
spect considering ’’granular structure" in the combustor is necessary. Existing evi- 
dence (ref. 2, fig. 8) indicates that the hot luminous portions of a rocket combustion sys- 
tem (as considered here) do not form a completely continuous medium, but that a consid- 
erable tendency exists for the overall flow system to be made up of pockets of hot gases. 
The time-dependent structure of such granulation depends on the jet breakup and mixing 
conditions prevailing in the combustor and is coupled with the burning itself. Conse- 
quently, conclusions regarding local conductivity (and possibly other properties) should 
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probably not be made on the basis of measurements for the overall volume of reacting 
gases. 

Finally, it is necessary to estimate a lower limit for the ion population produced in 
the system, which was previously used to determine an upper limit for the Debye shield- 
ing length. At the oxidant-fuel weight ratio and chamber pressure used, the following 
equilibrium composition is estimated (with 100 percent vaporization and burning efficiency 
assumed) (ref. 8): 


Species 

Mole fraction 

H 

0.019 

H 2 

.49 

h 2 o 

.47 

o 

.0002 

°2 

. 00008 

OH 

.006 


When each atom and molecule is assumed to act as a charge carrier, this composition 
gives 47 equivalents per pound. The present calculation is based on a view of electric 
current flow that depicts charges being carried in the flow stream by association with 
molecules. This problem involves overall mass transport and is to be distinguished from 
current flow between electrodes in a sealed gas-filled tube. The latter is a diffusion 
problem that is usually treated in terms of a sum- over- charge carriers with specified 
mean collision times (ref. 4, ch. 7). For currents used here, 

( — - — fara day\ ( Q 06 tQ 0 36 c\ _ 6x10 -7 to faraday/lb 

\96 500 C /V lb) 

With 1 faraday — - 1 equivalent (with equivalent as used here), 6x10”* to 4xl0”° equiva- 
lents per pound weight flow served as charge carriers, with only single ionization as- 
sumed. Dividing these values by 47 equivalents per pound gives an "effective degree of 

o _ o 

ionization" during current flow of x = 1x10 to 8x10" . With the temperature range of 
burning taken as 3000° to 3500° K at a pressure of 14 atmospheres, the ideal gas law 
gives the concentration of charge carriers as 57. 0 moles per cubic meter at 3000° K to 
49. 0 moles per cubic meter at 3500° K. This concentration gives an ion population of 

3 2 S 3 

N q = (moles/m ) (6. 0x10 molecules/mole) (x)(l ion/molecule) = ions/m 


where 
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N Q = 2. 9X10 17 

(for 

= 4. 5X10 17 

(for 

= 2.4X10 18 

(for 

= 2.7X10 18 

(for 


T = 3500° K; x = lXlO -8 ) 

T = 3000° K; x = lXlO" 8 ) 
T = 3500° K; x = 8xl0~ 8 ) 
T = 3000° K; x = 8xl0' 8 ) 


When these ion populations are compared with table I (p. 8), they are seen to be close to 
thermal ion populations that would require a burning temperature of 4000° K. Since the 
actual c* efficiency was only around 50 percent, there were fewer equivalents per pound 
than calculated here. Therefore these values for N q are probably low and represent 
only lower limits for the ion population in the hot portions of the combustion system (which 
are the parts to control). Consequently, the Debye shielding lengths calculated from 
these values represent maximum lengths for the system considered. 
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APPENDIX B 


SYMBOLS 


B magnetic induction (mean 

magnetic field in matter) 

C Coulomb 

c* characteristic exhaust velocity 

E electric field 

H Hartmann number 

h Debye shielding length 

I radiation intensity 

j vector electric current density 

k Boltzmann’s constant 

L characteristic length of Q 

1 vector element of length on a 

fluid line 

M general designation for atom 

or molecule 

N magnetic interaction parameter 

N q ion population 

o/f oxidant-fuel weight ratio 

P pressure 

Q volume element of fluid 


q electric charge 

R universal gas constant, (82.057 

(cm'*) (atm)/ (°K) (mole)) 

Re Reynolds number 

Re m magnetic Reynolds number 

T temperature, °K 

t time 

t characteristic electrical relaxa- 

c 

tion time in a fluid 
u vector velocity 

v^ ionization potential 

x degree of ionization 

€ permittivity of free space 

p magnetic permeability 

77 Avagadro’s number 

v kinematic viscosity 

p fluid density 

ct electrical conductivity 

u> frequency defined as 1 /t 
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